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Abstract. We search for signatures of high-frequency oscillations in the upper solar photosphere and low chromosphere in 
the context of acoustic heating of outer stellar atmospheres. We use ultraviolet image sequences of a quiet center-disk area 
from the Transition Region and Coronal Explorer (TRACE) mission which were taken with strict cadence regularity. The latter 
permits more reliable high-frequency diagnosis than in earlier work. Spatial Fourier power maps, spatially averaged coherence 
and phase-difference spectra, and spatio-temporal (k h ,f) decompositions all contain high-frequency features that at first sight 
seem of considerable intrinsic interest but actually are more likely to represent artifacts of different nature. Spatially averaged 
phase difference measurement provides the most sensitive diagnostic and indicates the presence of acoustic modulation up to 
/ « 20 mHz (periods down to 50 seconds) in internetwork areas. 
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1. Introduction 

In this paper we continue studies of solar atmosphere oscilla- 
tions based on analyzing the temporal brightness modulation in 
image sequences taken with the Transition Region and Coronal 
Explorer (TRACE) in ultraviolet passbands which sample the 
upper solar photosphere and low solar chromosphere. We again 
exploit the absence of seeing in TRACE data (apart from space- 
weather particle hits) to provide extensive Fourier diagnostics 
for quiet-sun network and internetwork areas with excellent 
sampling statistics. 

these tech- 



In iKriigeret all d2001 



henceforth iPaper IT) , 
niques were used in a comprehensive overview of quiet-sun 
brightness oscillation properties derived from TRACE im- 
age sequences in its three ultraviolet passbands centered at 
A = 1700, 1600, and 1550 A. In standard models of the so- 



lar atmosphere such as FALC of iFontenla et al.l i 1993b these 
passbands sample layers just below, at, and just above the 
temperature minimum, r espectively. The subsequent paper by 
Rutten & Krii gen ([2003) compared low-frequency ultraviolet 
brightness modulation at these wavelen gths to th e under lying 
white-light patterns in quiet-sun areas. iMuglachl 020031) ana- 
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lyzed TRACE ultraviolet brightness modulation maps contain- 
ing an active region. 

In this paper we return to the high-frequency a spects of ul- 
traviolet brightness modulation. The data used in IPaper II suf- 
fered from irregular timing intervals between successive im- 
ages, severely reducing the high-frequency information con- 
tent. The sequences used here have strict sampling regularity 
and are therefore better suited to search for high-frequency os- 
cillation signatures. We also employ a much improved align- 
ment method. 

The obvious motivation for such searches is given by the 
long quest for acous t ic hea ting of outer stellar atmosp heres 
started by iBiermannl dl948h and ISchwarzschildl dl948l) . It is 
concisely summarized by IWunne nberg et al. (2002), to whom 
we refer for further background. IWunnenberg et al.l 020021) em- 
ployed image sequences from the Gottingen Fabry-Perot spec- 
trometer at the German Vacuum Tower Telescope on Tenerife, 
scanning the non-magnetic Fe 1 5434 A line which samples lay- 
ers around h = 500 km above the white-light surface. They 
inferred the presence of sufficient power with 50-100 s peri- 
odicity (10-20 mHz in frequency) to compensate the radia- 
tive losses of the chromosphere, with apparent spatial power 
concentration above intergranular lanes. In this analysis we use 
TRACE data to search for corroborative evidence in ultraviolet 
brightness modulation from the same layers. 
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Fig. 1. Corrections for residual image displacements for the 
1600- A sequence, in the solar X (upper panels) and Y direc- 
tions (lower panels) 



2. Observations and data reduction 

plotted against frame number. The enlargements in the right- 
hand panels are for the short segments specified by the bars at 
left. 

The TRACE mission is described by Handy et al. (1999). 
We use ultraviolet image sequences downloaded from the 
TRACE archiveOl They were recorded on June 1, 2003 at the 
request of M. Carlsson (Oslo), who suggested strict cadence 
regularity and low data compression in order to minimize high- 
frequency artifacts, in particular those arising from timing ir- 
regularities as analyzed in Sect. 5 of IPaper J . TRACE was pro- 
grammed to obtain such image sequences in its 1600- A and 
1700-A ultraviolet passbands from 8:14 to 18:34 UT. We se- 
lected uninterrupted subsequences of 1120 images starting at 
11:23:12 UT and ending at 15:25:56 UT. They have strictly 
regular cadence at 13 s sampling interval in both passbands. 
The corresponding Nyquist frequency is fy y = 38.46 mHz; the 
frequency resolution is Af = 34.34 /iHz. The images sample a 
quiet area of 256" x 192" centered at X = -2.78", Y = 13.46" 
near the center of the solar disk, corresponding to a field of 
view of 512 x 384 square 0.5" pixels. The 1600-A and 1700- 
A images were alternately exposed for 1.724 and 4.872 s, re- 
spectively. The mid-exposure delay between the closest pairs 
of 1600-A and 1700-A images is 5.574 s. 

The image sequences were processed with the SolarSoft 
routine trace_prep described in the TRACE Analysis GuideB 
It corrects missing pixels (of which there were none in these 
data), replaces saturated pixels with values above 4095, sub- 
tracts the dark field, and corrects for the flat field. The dark and 
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Fig. 2. Brightness histogram of one 80-minute average of the 
1600-A sequence. The dotted lines define the split between net- 
work (right), intermediate (middle), and internetwork (left). 



flat fields were recorded on November 8, 2001 and March 11, 
2003, respectively. The image brightness was normalized by 
the exposure time. 

In the course of this analysis it became clear that pre- 
cise image alignment, including corrections for differential 
solar rotation and for spacecraft pointing jitter, is crucial to 
Fourier phase-difference analysis at high frequencies, and that 
we should considerably improve on the method used in lPaper I . 
In that paper, co-aligned subfields of the 1700-A sequences 
were cross-aligned to the corresponding subfields in the 1600- 
A sequences. This procedure copies alignment errors from 
one sequence to the other and so introduces a high-frequency 
phase-difference signal at the retardation set by the timing off- 
set betwe en the ex posures at the two wavelengths. It emerges in 
Fig. 1 8 of Paper J as a high-frequency drift of the spatially aver- 
aged phase-difference curves towards the average offset caused 
by non-si multane ous sampling shown in the center panel of 
Fig. 28 of lPaper II T he slowe r cadence of the October, 14 1998 
data also analyzed in lPaper II caused correspondingly larger off- 
set (Fig. 26). 

In the present analysis such erroneous cross-alignment sig- 
nals are reduced by significantly improving the alignment pro- 
cedure. In order to minimize the use of interpolation, we mea- 
sured pointing displacements per image through an elaborate 
procedure detailed below and then used these displacements to 
resample the original images directly onto an aligned grid. 

We began by shifting every row of each image in solar 
X to correct solar rotat ion including its dif ferential shear, us- 
ing the expression of lHowardetaiTd 19901) . We then aligned 
each image of 40-image 1600-A sub-sequences to the last one 
of the previous set, comparably to the procedure in Paper J . 
Each 1700-A image was subsequently cross-aligned to the cor- 
responding coarsely aligned 1600-A image taken 5.574 s be- 
fore. We then applied spatial smoothing through 5x5 pixel 
boxcar averaging to every image, and, merging the two se- 
quences, applied temporal smoothing per pixel by an eighteen- 
image boxcar average. Alignment of each individual image of 
the de-rotated sequences to this smoothed average yields dis- 
placement vectors per image. SolarSoft routine tr_get_disp 
was employed in all alignment computations. 
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Fig. 3. First panel: sample image from the 1600-A sequence taken at 11:32:22 UT. The intensity was clipped and scaled loga- 
rithmically in order to gain contrast in the internetwork. Second panel: 80-minute 1600-A average (12:44- 14:05 UT) using the 
same gray scale. Third panel: the pixel masks applied to all images in Sect. 13.21 Dark gray, light gray, and white respectively 
denote internetwork (28413 pixels), intermediate (36914 pixels), and network (7430 pixels). Black pixels are discarded. The box 
specifies the subfield selected for the power maps in Fig. [5] 



Figure Q] shows the resulting displacement corrections for 
the 1600-A sequence. These are the frame-by-frame residuals 
after the initial correction for differential rotation. They primar- 
ily describe pointing errors. Both X and Y components show 
oscillatory behavior with about 1.5-pixel amplitude and ap- 
proximately 100-minute periodicity caused by the spacecraft's 
orbital motion. The enlargements at right show ragged excur- 
sions with quarter-pixel amplitudes on short timescales which 
reflect pointing jitter. Solar rotation causes a much larger addi- 
tional drift in the horizontal direction. The de-rotation correc- 
tion ranges from -73.33 pixels at the equator to -72.83 pixels 
at the bottom of the field of view. Our use of whole-field align- 
ment automatically corrects for any departures from the ini- 
tially applied rotation law except for those in differential shear. 
The error estimates of lHoward et al.l (11990) suggest that the re- 
maining shear errors are within 0.004 pixel over our range in 
solar Y and time. 

In the final step of the ali gnment procedur e, the sophisti- 
cated algorithm described by iDeForesj d2004l) is used to re- 
sample the original images onto a 432x384 pixel grid corrected 
for differential solar rotation, for spacecraft orbital motion and 
jitter, and for the re-mapping from planar to spherical coordi- 
nates. The area of incomplete sampling due to solar rotation 
is discarded, as is a vertical strip at the left-hand edge of the 
field of view which erroneously appears bright in one 1700- A 
image. The resulting images consist of 432 x 384 0.348 Mm- 
square pixels. 

For part of our analysis, i.e., the spatially averaged Fourier 
spectra p resented in Fig. [7] in Sect. 13.21 we followed the pro- 
cedure of Paper I to divide the field of view into internetwork, 
intermediate, and network areas through classification of the 
time-averaged 1600-A brightness per pixel. Temporal averag- 
ing increases the contrast between the rapidly changing inter- 
network brightness and the more stable network emission. The 
1600-A sequence was split into three parts of approximately 



80 minutes duration. Figure |2]displays the brightness distribu- 
tion after averaging over one 80-minute part. It has a Gaussian 
peak and an extended high-brightness tail. A pixel is classified 
as internetwork if in all three 80-minute averages its bright- 
ness remains below the left-hand dotted line, which is chosen 
near the three peaks. Pixels with brightness above the right- 
hand dotted cutoff in all three averages are classified as net- 
work. Pixels that fall between the two lines in all three averages 
are classified as intermediate. Pixels that change category be- 
tween averages are discarded. This category amounts to 56% 
of all pixels due to the long sequence duration. It is large to 
avoid any mixing of internetwork, intermediate, and network 
behaviour. 

Particle hits were not corrected by interpolation but were 
removed on the basis of their high-frequency signature for the 
analysis in Sect. 13.21 Their single-image appearance produces 
anomalously strong high-frequency power. We applied a spatial 
mask to remove all pixels as well as their immediate neighbors 
that show Fourier power in excess of three times the average in 
the highest 50 frequency bins (36.7-38.5 mHz). 

Figure [3]presents a sample 1600-A image, one of the three 
80-minute averages, and the three masks. 



3. Analysis and results 

Fourier power, coherence and cross-power spectra were com- 
puted per pixel in the 1600-A and 1700- A image sequences 
over thei r full 24 3 -minute duration following the recipes in 
Sect. 3 of IPaper l We follow [Paper II also in the presentation of 
the resulting power maps, spatially-averaged temporal Fourier 
spectra, and two-dimensional (kh,f) diagrams. Here, the em- 
phasis is on high-frequency behavior and its significance. 
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3.1. Spatially resolved Fourier power maps 

As in IPaper il we distinguish three different normalization 
choices in displaying Fourier power per pixel as spatially re- 
solved maps, namely plotting the non-normalized oscillatory 
energy itself ("power"), 



P E (x,yJ) = \I(x,y,f)\ 2 



(1) 



the fractional modulation signal obtained by dividing the oscil- 
latory energy by the zero-frequency power ("modulation"), 



P f (x,y,f) = 



\l(x,y,f)\ 2 
\I(x,yM 2 



(2) 



and "Leahy" normalization obtained by dividing the energy by 
the zero-frequency amplitude, 



Pdx,y,f) 



\I(x,y,f)\ 2 
\I(x,y,0)\ 



(3) 



where x and y are spatial coordinates, / is the temporal fre- 
quency, and I(x, y, f) denotes the Fourier transform of the in- 
tensity measured by TRACE at location (x, y) at frequency /. 
Leahy normalization is used in the literature to estimate power- 



peak significance (e.g., Leah y et al 
but was not used in Paper it 

We here add 95% significance 
Grofhl d 1975b and 
of randomization 



1983, Dovle et al. 1999) 



estimation following 
first compare this to Fish e r's m ethod 
described b y, e.g., iBradlevI (11968b and 



Linnell Nemec & Nemec ( 1985) and applied to solar data by, 



appl: 

e.g. JO'Shea et al.ld200lb and lMuglacht (|2""003I) . Its assumption 
is that there is no signal at any frequency, so that the tem- 
poral order in which the data were taken becomes irrelevant. 
Comparison of the actual power spectrum to the spectra of 
temporal permutations of the data sequence then yields a sig- 
nificance estimate. The test is used iteratively, progressively 
deleting significant peaks until no new peaks are found. It puts 
no constraint on the noise power distribution at any given fre- 
quency, but the assumption that all samples are temporally 
uncorrected implies frequency-independent white noise. For 
large data sets it becomes impractical to repeatedly compute all 
possible permutations. Actual tests are therefore usually lim- 
ited to a few hundred permutations, but even then remain com- 
putationally expensive 



The much simpler significance estimation of Groth ( 1975) 
assumes that at any frequency the real and imaginary parts of 
the Fourier power have independent normal distributions. It re- 
quires explicit specification of the noise power as a function of 
frequency, i.e., the noise is not assumed to be white. 

We compare the randomization test with Groth's test in 
Fig. E]for an internetwork and a network pixel, adopting 95% 
confidence levels in both tests. In the randomization test, a 
power peak that is above the maximum power of the random- 
ized data in more than 95% of 500 permutations is consid- 
ered statistically significant. Such peaks are subsequently re- 
moved in the iterative re-application of this procedure until no 
more significant peaks are found. For Groth's test we decided 
from visual inspection to assume that the power spectra dis- 
play white noise above / = 24 mHz. The corresponding 95% 
significance cutoff is 2.996 times higher than this noise level. 
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0.5 0.45 




4 8 12 16 20 24 28 32 36 
frequency [mHz] 

Fig. 4. Comparison of the random ization test to significance es- 
timation following iGrothl (1 1975b . Upper panel: network pixel. 
Lower panel: internetwork pixel. The ragged curves are the 
temporal Fourier power at 1600 A using TRACE data units 
divided by the exposure time, on logarithmic scales. In each 
panel, the top row of tick marks identifies all significant peaks 
according to the randomization test. The second and third rows 
identify significant peaks at the 99.9999% and 95% signifi- 
cance levels using Groth's test assuming white noise and ab- 
sence of signal above / — 24 mHz. The dotted lines show the 
corresponding cutoff levels. 



The top and bottom rows of tick marks in Fig. |4] specify 
the positions of the peaks that are estimated to be significant 
by the two methods. It is obvious that the randomizat ion test is 
far mo re rigorous than Groth's test, as pointed out bv lMuglachl 
(2003). The middle rows of ticks result when the upper dot- 
ted line is used as Groth-test cutoff level, at fourteen times the 
noise level corresponding to 99.9999% confidence. It closely 
matches the peak-finding by the randomization test. Thus, for 
these data a stringent Groth test may replace the randomization 
test at much smaller computational cost. This is likely to hold 
for other data with white noise. 

The network pixel in the upper panel of Fig. |4]has larger 
low-frequency signal than the internetwork pixel in the lower 
panel, a different power hump around 5-min periodicity, and 
higher high-frequency noise but rather similar peak survival 
above the lenient 95% Groth cutoff estimate. 

Figure [5] expands such comparison of internetwork versus 
network by displaying spatial power maps for the small but il- 
lustrative subfield specified by the rectangle in Fig. [3] for both 
the 1600-A and 1700-A passbands. The subfield is shown in 
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Fig. 5. Spatially resolved power maps using different methods of normalization for the subfield shown in Fig. [3] The grayscale 
displays the logarithm of the temporal Fourier power, clipped to improve contrast. Columns: different frequency bands as specified 
at the top. Rows: 1600-A and 1700-A passbands, with different power normalization as specified in the first column. 



different temporal frequency bands, with the three different 
normalizations (power, Leahy, and modulation, respectively), 
and finally also without normalization but with all pixels hav- 
ing power below the 95% Groth cutoff made black. 

The first column shows low-frequency power. It shows the 
stable nature of bright network. The power-normalized modu- 
lation maps show noisy behavior from the division because the 
frequency band is close to zero Hz. 

The second-column frequency range of 2.6-3.6 mHz cor- 
responds to periodicities around 5 minutes. The network ap- 
pears power-bright in the unnormalized maps, about equal 
to the internetwork in the amplitude-normalized Leahy maps, 
and power-dark in the modulation maps. Thus, the choice of 
normalization affects the apparent relative dominance of net- 
work an d internetwork oscillations, as discussed extensively in 
Paper l Note that in all representations a power-dark moat ap- 
pears around the network. 

The 5-7 mHz maps describe the chromospheric three- 
minute oscil lation which pervades internetwork areas (e.g., 
Rutten 1995b . They indeed show the network power-dark in 



all representations. There are irre gular power-bright "aureole" 
patches near network (cf. Pap er J) . 

The two high-frequency columns on which we concentrate 
here illustrate the care that must be taken in interpreting such 
power maps. In the rightmost column (28-32 mHz) the net- 
work stands out very brightly in the unnormalized power maps, 
inviting a claim that the magnetic elements making up the net- 
work display high-frequency wave heating. On the other hand, 
the network appears power-dark in the modulation maps, invit- 
ing a claim that high-frequency oscillations are suppressed in 
magnetic elements. However, the close spatial correspondence 
of both these bright and dark features with the bright network 
in the unnormalized maps in the first column suggests strongly 
that they are simply due to the larger overall network bright- 
ness. A similar power-contrast flip is seen in the 12- 16 mHz 
maps for the 1700-A passband, but the unnormalized 1600-A 
map for these frequencies appears rather featureless. The latter 
copies directly into the Groth map, but with considerable pixel 
deletion wherever the power averaged over the frequency range 
drops below the cutoff. In the rightmost column the Groth maps 
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accept only a minor fraction of the pixels as significant, both for 
the network and the internetwork. A more strict criterion, and 
certainly the 95% iterative randomization test, would reject all. 
The patterns seen in the other 28 - 32 mHz maps thus are most 
likely artifacts caused by sources of high-frequency errors with 
some sensitivity to the low-frequency power. Note that Leahy 
normalization diminishes the apparent structure for 1700 A but 
turns it power-dark at 1600 A. 

The bright specks in the 28 - 32 mHz maps are due to par- 
ticle hits. They produce high-frequency signal through their in- 
stantaneous appearance. The pattern of horizontal stripes re- 
sults from TRACE'S JPEG data compression (K. Muglach, pri- 
vate communication). It appears as a grid pattern with 8-pixel 
mesh size in comparable power maps computed from the orig- 
inal non-aligned image sequences. The compensation for solar 
rotation smears out the vertical grid components, leaving only 
the horizontal ones. 



3.2. Spatially averaged Fourier power, phase 
difference and coherence 

We now turn to temporal Fourier analysis with spatial averag- 
ing over the different pixel categories defined by the third panel 
of Fig. [3] The averaging is performed on the Fourier measure- 
ments per pixel. It reduces the noise in these measurements and 
therefore improves the dete ction of relatively small modulation 
signals. Following Paper I, the power and coherence are aver- 
aged directly over all relevant pixels. In the present analysis we 
compute coherence per pixel using frequency smoothing over 
9 bins rather than 5. 

The phase differences ar e again averaged with cro ss-power 
weighting as introduced by iLites & Chipmanl dl979l) . i.e., the 
spatial average over the phase differences of all pixels trans- 
mitted by the mask per frequency bin is defined as the angle 
of the vector sum of the cross-powers of all contributing pix- 
els with reference to the real axis. The advantage of such vector 
averaging is that it makes signals stand out even in the presence 
of much larger noise. For pure noise the vector mean does not 
go to zero or some other definite value but fluctuates randomly 
over the full -180 to +180 degree range between adjacent fre- 
quency bins. A small signal, much smaller than the noise, may 
therefore emerge as a systematic pattern across multiple bins. 

This is illustrated by Fig. [6] which displays simulation 
results for vector-averaged cross-power distributions of pure 
noise (first panel) and of pure noise with a much smaller su- 
perimposed signal with fixed phase difference (second panel). 
In the latter case, the vector addition over 30000 pixels with 
small but systematic signal shifts the much wider scattercloud 
to a location well separated from the origin, making the phase- 
difference distribution in the third panel strongly peaked. 

Figure |7]presents results from the new TRACE sequences 
in the form of power, coherence, and phase-difference spec- 
tra with spatial averaging separated between the internetwork, 
intermediate, a nd ne twork areas. The format is similar to 
Figs. 18-22 of IPaper II but adds two dotted lines in the phase- 
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Fig. 6. Illustration of vector-averaging phase differences for 
30000 pixels. First panel: distribution of the cross-power vec- 
tor sum for pure Gaussian noise in the complex plane. Second 
panel: same as the first panel, but with a signal with amplitude 
of only 3% of the rms noise with degree phase difference 
added. The vector summation of the 30 000 samples shifts the 
scattercloud significantly to the right. Third panel: correspond- 
ing phase-difference distributions for pure noise (dashed, flat) 
and with the signal added (solid, peaked). The latter reaches 
18.4 at phase difference with 20-degree half- width. 



sures in the two passbands. We have aligned the two sequences 
to match the closest pair combinations, corresponding to the 
upper dotted lines. All of our phase-difference evaluations are 
corrected for this sampling offset, which means that signals 
that are intrinsically in phase in the two passbands should in- 
deed end up along the horizontal A(f> — axis. Measurements 
that end up on a dotted line imply modulation with phase de- 
lay exactly matching the corresponding sampling delay. The 
grayscaled scattercloud represents individual pixels. In the case 
of pure noise the lcr rms estimates cover 68% of the full figure 
height around a randomly fluctuating mean. 

The results in Fig. [7] are similar to those in Fig. 18 of 



difference panels. These are the phase shifts associated with 
the temporal sampling delay due to the non-simultaneous expo- 



Paper I except for the high-frequency phase differences of in- 
terest here. The present results are more reliable thanks to the 
regular sampling cadences, lower data compression, and better 
image alignment. 

The internetwork phase differences reach a wide maximum 
at / » 7 mHz and then remain well defined at positive values 
up to the Nyquist frequency, but with increasing noi se above 
20 mHz. There is no drop to negative values as in IPaper I , 
which we now attribute to the cross-alignment used there as 
discussed in Sect. [2] above. However, the present results also 
show a drift to the phase difference associated with the tim- 
ing delay at high frequencies. The internetwork power spectra 
show acoustic humps around 4 minutes and become negligible 
above 12 mHz. The coherence also peaks around 4 mHz and 
drops to the 9-bin noise level near 20 mHz. 

The network phase differences in the third column of Fig. [7] 
are much noisier due to the far smaller number of pixels. 
Nevertheless, they show a narrow peak of increased phase dif- 
ference and reduced coherence aro und thre e-minute periodicity 
which is not present in Fig. 18 of Paper J and which we deem 
significant. At high frequencies they become more erratic and 
shift to the timing correction line, which points to a systematic 
error. 
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Fig. 7. Temporal Fourier spectra, spatially averaged over internetwork (left), intermediate (middle), and network (right), plotted 
against frequen cy up to the Nyquist limit with the corresponding periodicities shown along the top. The format corresponds to 
Figs. 18-19 of Paper J . adding dotted lines indicating the sampling time delays and omitting lcr rms estimates for power and 
coherency to avoid clutter. Upper panels: phase-difference spectra. Lower panels: coherence (upper curves) and power spectra 
(solid for 1600 A, dashed for 1700 A). The random-noise estimate for the coherence is C — 0.33. The power spectra are on linear 
scales and are all scaled by the same factor. 



The intermediate-class pixels in the center column of Fig. Q 
produce primarily internetwork-like behavior. 

3.3. Two-dimensional Fourier power and phase 
difference 

Figure [8] presents two-dimensional Fourier power and phase 
difference spectra in the form of (k^,/) diagrams. They mix 
the network, internetwork and intermediate areas. Particle hits 
were not removed because their contribution to the noise is 
small except at high spatial and temporal frequencies where 
the diagrams are noisy anyhow. 

The power and phase differences are averaged over rings 



of constant k^. 



with k\ 



k\ + ki, assuming absence of pre- 



ferred horizontal propagation directions. The number of sam- 
ples per ring increases with k^ up to the Nyquist frequency 



per axis k x Ny 



= k 



y.Ny 



= 9.0 Mm . Beyond this value, k^ 



can still be computed but with fewer samples and increasing 
loss of isotropy in each successive bin, up to k^ = V2£ X N y = 
12.8 MirT 1 which samples only oblique propagation. 

The left-hand panel of Fig. [8] shows the (k^, /) diagram for 
1600- A power. The acoustic p-mode ridges and pseud o-ridges 
above the Lamb line were extensively discussed in iPaper I . 
There is no particular structure evident in the high-frequency 
regime of interest here. At low frequencies there is a ridge of 



enhanced power at high spatial wavenumbers, approximately 
corresponding to / = (1/2^)2 kms -1 which is caused by the 
compensation for solar rotation. Features that are fixed to the 
CCD camera, such as the results of an imperfect flat field, or 
"hot" pixels, move apparently with this speed against the di- 
rection of solar rotation and produce power. 

The right-hand panel of Fig. [8] displays the corresponding 
(fch,/) diagram for phase difference between the 1700- A and 
1600- A sequences. The acoustic ri dges sta nd out through larger 
phase difference, as discussed in IPaper j The wedge of nega- 
tive phase difference at low frequencies and wave numbers was 
attributed to atmospheric gravity waves by iRutten & Kr ijger 
(2003). The effects of solar rotation are again visible as a ridge 
of slightly increased phase difference. We attribute this increase 
to the systematic (-0.285, 0.273)-pixel image offset between 
the two passbands before alignment, which causes apparently 
traveling features fixed to the CCD to appear at a given solar 
location with some time delay in the two passbands. 

At the highest temporal frequencies of interest here, noise is 
easily identified as pepper-and-salt patterning where the phase 
differences jump widely from one bin to the next. A large patch 
of smooth variation extends up to about 20 mHz in frequency 
and 4 Mm 4 in wavenumber. This patch contributes most to the 
definite phase behavior in Fig. [7] which is therefore set by these 
spatial scales. 
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Fig. S.Left: power in the 1600- A image se- 
quence plotted as function of horizontal 
wave number and temporal frequency 
/. The logarithmic grayscale is clipped to 
show the ridges around 1 Mm" 1 and 5 mHz. 
The slanted line is the Lamb line / = 
(l/2n) c s kh with c s = 7 kms~'. Right: cor- 
responding phase differences between the 
1700-A and 1600-A image sequences. To 
avoid figure clutter, contours are only shown 
if they lie below the dashed curve or en- 
close a large area. The grayscale is clipped 
at -9 and 45 degrees to increase contrast. 
The white blob peaking at k h = 9 MirT 1 and 
/ = 5 mHz reaches 110-degree difference. 
The pepper-and-salt regions reflect noise. 



The conspicuous white blob of large positive phase differ- 
ence located at 0.75 Mm wavelength and 3-minute periodicity, 
with an extended tail upward, is enigmatic. The left-hand dia- 
gram suggests enhanced power at this location. It seems likely 
that the peak in the network panel of Fig. |7]corresponds to this 
blob, and that therefore the source should be sought in the net- 
work. It is very tempting to attribute it to solar three-minute 
waves with small horizontal extent in magnetic elements. Its 
exceedingly large value (up to 1 10 degree difference at its cen- 
ter), the clear reduction in coherence, and the lack of such 
a blob in comparable phase-difference diagrams sampling the 
Ca ii H line core and inner wing from observations made with 
the Dutc h Open Telescope (unpublished analysis analogous to 
Fig. 7 of iRutten et al . 2004), would then suggest phenomena in 
the transition region that leave a signature in these data through 
the C iv doublet at A = 1548 A and 1550 A (cf. lHandv et al ' 



1998). However, the blob lies at the spatial Nyquist frequency 
per horizontal axis, and its shape varies with changes in the im- 
age alignment procedure. We reluctantly conclude that the blob 
is likely a TRACE artifact, presumably of instrumental origin, 
introduced by the data processing, or a combination of these. 

4. Discussion 

We find intriguing high-frequency behavior in all our Fourier 
displays, but at the same time also find reasons to disbelieve 
these signatures above 20 mHz or even lower frequencies. 
The pixel-by-pixel power maps in Fig. [5] have disconcerting 
high-frequency contrast sensitivity to the type of normaliza- 
tion above 10 mHz. The spatial averaging in Figs. [7] and [8] re- 
spectively over mask-selected pixel types and annuli, improves 
the significance in phase-difference measurement, but the high- 
frequency behavior in Fig. [7] is puzzling in the trends towards 
the timing delay lines and the absence of purely random behav- 
ior even at the highest frequencies. The prominent white blob 
in the phase-difference panel of Fig.[8]is presumably an artifact. 



The phase-difference averaging with cross-power weight- 
ing over the different pixel categories employed in Sect. 13.21 
is the most sensitive method to identify weak oscillation sig- 
natures in the presence of noise. Each phase-difference dia- 
gram in Fig. [7] indicates systematic non-random behavior out 
to frequencies far beyond the extent of measurable power or 
even of measur able coh erence. This was already the case in 
Figs. 18-19 of Paper II and also in the similar diagra ms from 
groundbased Can H spectrometry in Figs. 20-22 of iPaper I . 
It is attractive to believe that the cross-power weighting indeed 
enhances the sensitivity of the phase-difference measurement 
to very small signals otherwise drowned in noise out to well 
above the coherence limit of at most 20 mHz, but it is alarming 
that even at the highest frequencies the phase differences do not 
show the randomness expected, and that in all three panels in 
Fig.|7]it seems to favor the instrumental timing correction. It is 
likely that residual image alignment errors are the cause of this 
anomalous behavior. 

It is well-known that the increasi ng lack of response due 



to wide c ontribution functions (e.g ., Beckers & Milkev 11975 
Deubner|[l976[ ISchmieder & Meinlll980l |Purranllll98ri ham- 



pers the detection of high-freque ncy signals. It was recentl y 
elaborated in TRACE context by IFossum & Carlssonl (120041) . 
In addition, we have learned from M. Carlsson (private com- 
munication) that simulations of acoustic waves propagating up- 
ward in the solar atmosphere as in the we ll-known Can H2V 
grain simulation of ICarlsson & Steinl (1 19971) subjected to com- 
putational 1700- 1600-A phase-difference analysis which em- 
ulates the observational analysis presented here meets unex- 
pected computational problems at low signal to noise and high 
frequencies. 

On the other hand, we have reproduced our phase- 
difference results in tests using double precision computation. 
Very similar non-random positive phas e-differe nce behavior 
also appears up to 20 mHz i n Fig. 20 of Paper I, based on the 
Can H slit spectrometry of ILites et al. (1 19931) and measured 
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from Ca n H wing intensities and Fe i blend Doppler shifts 
formed at lower and similar heights as the ultraviolet continua 
sampled by TRACE. The comparable signature in /-/ and V-V 
diagnostics with negative V-I lag shown there is in agreement 
with acoustic waves. The steep V-V signature of upward prop- 
agation in internet work areas present in the lower-left panel of 
Fig. 21 of IPaper II seems significant also. 

In summary, the coherence spectra in Fig. [7] the close 
agreement of the phase diffe rences i n that figure with those 
from Can H spectrometry in Paper it and the smoothness of 
the corresponding gray area in Fig. [8] all taken together lead us 
to believe that the phase-difference signals derived from these 
new TRACE sequences have a solar origin up to 20 mHz at 
least in the internetwork, and are to be attributed to acoustic 
waves. 

This conclusion supports the detec ti on of high-frequency 
acoustic waves by IWunnenberg et al .] d2002h as significant 
Doppler-shift power in the 10-20 mHz frequency band from 
differential Fe i measurements addressing similar atmospheric 
heights. The drop of power with frequency in our Figs. [4] [7] and 
[8] suggests that their detection is dominated by the lower fre- 
quencies in this band. Our results indicate wave presence also 
at the higher frequencies. 

The ultraviolet continua used here suffer from strong scat- 
tering while the TRACE filter bandwidths are wide and overlap 
considerably. Numerical simulations as those presently under- 
way at Oslo may explain how and why the phase differences in 
Figs.|7]and[8]level out at positive values. In our opinion, com- 
parison with detailed numerical simulations is also required to 
substantiate any claim that acoustic waves in the 10-20 mHz 
regime compensate fully for the radiative losses of the chromo- 
sphere. 

5. Conclusion 



New ultraviolet image sequences from TRACE give evidence 
of brightness modulation up to 20 mHz in quiet-sun internet- 
work. We interpret this signal as a signature of acoustic waves. 
It is similarly present in Ca n H and Fe 1 I— I and V-V phase- 
difference spectra in Fig. 20 of Paper II and it supports the detec- 
tion of acoustic wave power in the 10-20 m Hz frequency band 
from F e i Doppler-shift measurements by IWunnenberg et al. 
(2002). The evidence for modulation at higher frequencies re- 
mains inconclusive. 

TRACE-like ultraviolet imaging will be achieved with the 
Atmospheric Imaging Assembly on NASA's Solar Dynamics 
Observatory, but it is not yet clear whether its hardware and op- 
eration will permit better high-frequency modulation measure- 
ment than with TRACE. New ground-based telescope technol- 
ogy, in particular large aperture combined with adaptive optics, 
will provide accurate Doppler shifts from integral field spec- 
troscopy at high cadence and low noise of the same layers using 
appropriate spectral lines in the optical. Numerical simulations 
may contribute quantification of the corresponding energy bud- 
gets. 
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